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We study the one-loop radiative corrections for massless fermions in de Sitter space induced 
by a Yukawa coupling to a light, nearly minimally coupled scalar field. We show that the 
fermions acquire a mass. Next we construct the corresponding (nonlocal) effective fermionic 
action, which - in contrast to the case of a massive Dirac fermion - preserves chirality. 
Nevertheless, the resulting fermion dynamics is precisely that of a Dirac fermion with a 
■ mass proportional to the expansion rate. Our finding supports the view that an observer or 

a test particle responds to a scalar field in inflation by shifting its energy rather than seeing 
a thermal bath. 
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The question of the dynamics of fermions interacting through a Yukawa coupling with a massless 



O 
O 

(N 



^jT)| minimally coupled scalar field in de Sitter space was originally considered in Ref. lj. It is the 

simplest nontrivial example allowing for a one loop study of quantum effects during inflation. Here 
we reconsider the problem, by assuming that the scalar field is light and nearly minimally coupled, 
which allows us to calculate radiative corrections employing a de Sitter invariant scalar propagator. 
We make substantial analytic progress, such that some of the questions that were left unanswered in 
the original work, e.g. the problem of Pauli blocking and vector current conservation, are resolved. 

Gaining analytical understanding of fermion dynamics during inflation may have important 
cosmological ramifications, in particular for baryogenesis and the origin of dark matter in the 
Universe. Closely related are effects in massless scalar quantum electrodynamics, which may be 
of relevance for cosmic magnetogenesis P, Q] , and of graviton induced self energy jj] in de Sitter 
space. 

On a more conceptual level, quantum theory in de Sitter space is most commonly considered 
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in the context of horizon thermodynamics, which is ocally witnessed by an Unruh detector at a 
response rate corresponding to thermal radiation j], f|. In recent work, we argued however that 
local effects predominantly become manifest in a shift of the detector's energy levels 3|. Since the 
Unruh detector is an idealised device, it is of interest what a test particle in a proper field-theoretic 
setting experiences, in particular whether effects similar to those of a thermal bath are present. 



2. ONE-LOOP SELF ENERGY 



Massless fermions are conformally invariant, such that in a D-dimensional spacetime with the 
metric 



ds 2 = a 2 (-drj 2 + cht 



(1) 



their propagator arises from a conformal rescaling of the flat-space counterpart 



iS(x; x) = (aa) 2 i<j 



47T 2 



Ax 2 (x; x) 



(#) (aa'^ij^AXf, 



2vrT [Ax 2 (x;x')] 



D_ 1 
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(2) 



where Ax 2 = —(\rj — rf\ — ie) 2 + ||x — x'|| 2 and Ax 
De Sitter space is endowed with the scale factor 



1 

llrj 



(3) 



where H denotes the Hubble parameter and r] conformal time. The Green function for a light 
nearly minimally coupled scalar in de Sitter space is of the following de Sitter invariant form 



>D-1 



iAfa ,^ = W + ^)r(^-^ gC , 2 



(te)Tr(f) 



where 



D - 1 



,,D-1 D-l D 
2-f ? i( — h v D , — v D , — , 1 



(4) 



m 2 + £R D 



(5) 



2 / //-' 

m is the mass of the scalar field, Rr> = D(D — 1)H 2 denotes the Ricci scalar curvature of de 
Sitter space, and £ is the coupling constant of the scalar field to curvature (£ = 1/6 corresponds 
to conformal coupling, £ = to minimal coupling). The function y = y(x;x') is related to the de 
Sitter invariant length, i = £(x; x'), as follows: 



4 sin 2 (-Hi) = aa'H 2 Ax 2 



Ax 2 



(\rj — rj I — ie) + ||x — x 



/ 1 1 2 



(6) 
(7) 
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where a = a(rj) = —l/(Hr]), a' = a(r/) = — 

If m 2 + £,Rd = 0, the Green function (jlj) is ill-defined. When attempting to calculate it 
by performing the momentum-sum over the modes, one sees that this problem corresponds to an 

□ n 

infrared divergence |9|, which has to be regulated. There exists no de Sitter invariant regulator |10], 
which implies that there is no de Sitter invariant Green function for a massless minimally coupled 



3, 



the use of a de Sitter invariant 



scalar field. In a work on scalar electrodynamics in de Sitter space 
propagator with a small mass or weak coupling £ to gravity was used in studying radiatively 
induced photon mass generation in de Sitter space. Here, we adapt this method to study the 
fermion dynamics in de Sitter background. 

Since we are primarily interested in considering a small scalar mass and weak coupling to gravity, 



m 2 + iR D ^H 2 , (8) 



we find the following representation of the scalar propagator useful ( cf. also Ref . 



n /„\ n 



(4vr)f r(i + ^)r(i-^)l^ r(f + n)r(n + i) 1 >! 



™ r(| + z^ D +n)r(|-i/ D +n 



n=-l 



r(3- § +n)r(n + 2) U 



-f+2- 



Notice that the conformal propagator corresponds to the n = — 1 term of the second series (con- 
taining the -D-dependent powers of y). We have recast the scalar propagator into the form @, 
such that the cancellation of all terms n > is manifest in D = 4. The cancellation occurs only to 
leading order in an expansion around D = 4 however. Since the prefactor is singular at D = 4, the 
true and finite result in D = 4 is obtained when the series are expanded to linear order in D — 4, 
resulting in coefficients that are suppressed by (in 2 +£Rd) I H 2 "C 1, as we show below. The explicit 
form for the propagator that is suitable for this problem is thus (y = [(3/2) 2 — (m 2 + £R)/H 2 ] 1 ^ 2 ), 

H D - 2 „sD \ 1 



-if)(l + nj - i){2 + nj + 0{D - 4) . (10) 
The fermionic and scalar propagators can now be assembled to give the self-energy 
- iS(x; x') = (-if^a D ) iS(x; x') (-i//x^a ,D ) iA(x; x') + i6Z 2 (aa') £ ^ 1 i@6 D (x - x') , (11) 
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where / is a dimensionless Yukawa coupling constant, 5Z 2 denotes the fermionic field-strength 
renormalisation counterterm and is a renormalisation scale. 
More explicitly, the one-loop fermion self-energy is given by 



S,T 



D 



Ax 



f 2 H 2 (aa')2 f 2 lN i i 7 ^Ax /j 

2 5 7T 4 

/ 2 tf 2 (aa') 



2D-2 

y 
4 



2 8 7T 4 



lw9 2 \i7 MAx 



Ax 2 



hi 



+ ... + i5Z 2 (aa')~i^ D (x - x') 



5 
" 2 

(12) 



where ip(z) = d[ln(T(z))]/dz, je = 0.57.. = —ip(l) is the Euler constant, 27^ — 1 = —ipiX) — i/}(2), 
2"fE — 5/2 = —ip(2) — and v = vd=a- Note that in (|12[). we kept the first term in general D 

dimensions, since this term requires renormalisation, while the other terms contain only integrable 
singularities, and hence in writing them we took the limit D = 4. The higher order terms we 
neglected in (fT2*)) are all suppressed as (m 2 + £ s R)/H 2 . 

Following Ref. [lj, we use the minimal subtraction scheme to extract the contribution from the 
first term in ((12)) leading to an ultraviolet divergence in D = 4. The counterterm is 

f r(f-i) 



5Z 2 



(13) 



(Z>-4)(£>-3)- 

Eventually taking D — > 4, we find the renormalised fermion self-energy (up to linear terms in 
(m 2 +ZR D )/H 2 ), 

f 2 (aa>) 



E(x; x') 



2 10 vr 4 
f 2 {aaf)t 



2 5 vr 2 
f 2 H 2 {aa')^ 



$ 8 4 [ln 2 (^ 2 Ax 2 ) - 2 ln(^ 2 Ax 2 ) 
■ln(aa')i^ 4 (x-x') 



+ 



2 9 



7 2 # 4 (aa') 



^Lilp - i) <9 2 ln 2 (iO/ 2 Ax 2 ) + 21n(aa')0 <9 2 ln(# 2 Ax 2 )} 



(V - -i/ 2 ) ln 2 (K!# 2 Ax 2 ) + 21n(aa')^ ln(tf 2 Ax 2 )} , (14) 



2 10 7T 4 



where we defined 



K 



1 3 3 

-exp 27s - 1 + V(- + v) + ^(g ~ ^) 

1 



4 6XP 



(15) 



When expanded around s = (m 2 + ^R)/H 2 = 0, we find 
A" 
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1 
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3H 2 1 


~ -exp 
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~ —exp 
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-cxp 



(16) 



Note that the second term in the self-energy Q14|) breaks de Sitter invariance anomalously. In 
the flat space limit, where H = and the scale factor a is constant, it can always be removed by 
either setting a = 1 or by an appropriate choice of the renormalisation bZ<i- In de Sitter space 
however, the momentum scale is continuously shifted, such that this term does not vanish at all 
times. 

3. EFFECTIVE EQUATION OF MOTION 

In order to study the dynamics of fermions, we shall need the retarted self energy, 



J ret 



(17) 



For the calculation of £ re t the following identities are useful (recall an additional minus sign in 

£+-), 



ln(aAx^_ H 



In 2 (a Ax 



2 ^ 
++J 



- ln(aAx + -)=2iir6(Ar] - r)6(Arf) , 
ln 2 (aAx + __)=4ivr6'(Ar/ - r)9{Arf) In \a(Ar] 2 



(18) 



where 6 = 9(x) denotes the Heaviside step function, 6(Ar] 2 — r 2 ) = 6(An — r), r = \\x — x'\\ and 
we used 



Ax 2 



(\Ar]\ - ie) 2 + \\x- x'\\ 2 , Ax 2 + _ = -(Arj + ie) 2 + \\x-x 



(19) 



Using this, we can now write the retarded one-loop fermion self-energy as 



''^-id d 4 \e(A v - r)0(An) [in |u 2 (Ar? 

7T I. L 
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/W) 

2 5 7T 2 

f 2 H 2 (aa' 



.2 r 2. 



ln(aa')i^ 4 (x - x') 



+ 



f 2 H A (aa'y. 

2 8 7T 3 



-) {i$d 2 [0(Ar] - r)O(Ari) In \KH 2 (Ar] 
+ ln(aa')# d 2 \0(Ari - r)9{Arj)] } 

-)(- 

4/V4 



2 r 2s 



[v 2 - t) (t - {tf [0(Ar/ - r)6{Arj) In \K x H 2 {Ar] 2 - r 2 



+ 



ln(aa')i# \o(Ar) - r)6(An)j } . 



(20) 



We shall now use the result (|20|) to solve the one-loop Dirac equation 1 

ak$JiP(x) - /dV£ ret (x;x>(z') = 0. 



Let us define the conformally rescaled wave function 



X (x) = a i/2 ^{x) = j k ' S X {ri) , 



(21) 



(22) 



which by spatial homogeneity allows for the decomposition 



x(v') 



(23) 



such that upon inserting the retarded self-energy ()20|) into (j21j) and performing the angular inte- 
grations, we get 

tyxfa) + -ir^idd 4 -e lk - x / drf / drrsm(kr) In l/^Ar/ 2 - r 2 )\ - 1 x {jf) 
2 b 7H k J Jo L J 

/ 2 



2 5 vr 2 



ln(a)i^x(x) + i$( ln(a)x(x) 

f 2 u~2 1 \ f 1 - f*7 

+ A__ (V _ j J $ d 2 - e ikS I drj' a! \ drr sm{kr) In [KH 2 (Arj 2 - r 2 )}x{r,') 



i - /-n /-At) 



+ ln(a)i$d-e J drfa! J drr s'm(kr)x(r]') 
+ i$d 2 -e ikx J dri'a'ln(a') J drr sm(kr)x(r]') 



fW (2 1 U 9 

2 6 7T 2 



1 - /"'7 /"At) 

' d?/o / (irrsin(A:r)ln[i^i J ff 2 (Ar ? 2 -r 2 )]x(V) 
o 



+ ln(a)i$ - e lk x J drj cl J drr sm(kr)x(rj') 
i- At] 



i$ J drj 1 a 1 ln(a') J drr sin(fcr)x(i/) 



(24) 



Upon performing the radial integration this becomes 



\0\(.r) / '(/,/<{ 21n(//A;/)- 1 sm(/,A//) - /, A//cos(/,A//) 



+(£A7 ? ) 2 £(£A7 ? ) x (r/) 



/ 2 



25tt 2 



ln(a)i^x(«) +i^(ln(a)x(^ 



/^a, 2 _1 
2 5 tt 2 V 4 



-!- -^ r- , [ ir - - j<ji^9 2 J_ e ifc-x / ' dr? ' a '<{ ln(KH 2 Ari 2 ) sin(/,A//j - AA,/cos(7,A//) 



1 Note that the sign in front of the nonlocal term differs from the sign in Ref. 0. The sign error in has been 
corrected in the corresponding erratum. 
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+(fcAr ? ) 2 ^Ar / ) \ x (rf) 



+ ln(a)i 



k 3 



fH^a 2 
2 6 



d?7 a sin(^Ar/) — kAr] cos(A;Ar/) 
drf a! \n.{al) s'm(kArj) — kAr] cos(kArj) x(?/)| 
(z^ 2 - i) (| _ ^ 2 ) JLjkzJ V d V 'a' 2 { In^H 2 Arf) [ sm(kAr]) - kAr] cos(kAr]) 



+{kAr,) 2 ttkAr]) \ x {rf) 



+ ln(a)i$ -g^e 



ik-x 



drj'a' 2 sin(kArj) — kArj cos(kAr]) xW) 



I drj'a' In(a') sin(kArj) - kAr] cos(kAr]) xW) } = , (25) 



where we made use of the following integrals, 

z 2 £(z)=z 2 / dxx sin(zx) ln(l — x 2 ) 
Jo 



(26) 



ir 7Tir i r / z \ 

2sin(z) — cos(z) + zsin(z) si(2z) H — + sin(z) — zcos(z) ci(2z)—je — In ^— J 



where 



si(z): 
ti(z)-- 



sin(i) 
t 

cos(t) 

~t 



dt 



dt 



sm(t) , 7T 

— —dt , 

o t 2 ' 



cos(t) — 1 

t 



dt + J E + M?) 



(27) 



Because the integrands vanish at the upper limit of integration at least as (A77) 3 ln(Ar/), the 
operator d 2 = —{d 2 + k 2 ) acting on the integrals in 1)251) commutes with J drj' , and may be directly 
applied to the integrands. The result is 



lh-X 



drj'\ 2ln({iAr/) sm(kAr/) — cos(kAr/) si(2kAr/) + ~ 



+ sin(kAr]) ci(2kArj) — je — In 



kAi] 



+ 



2 5 vr 2 



h\(a)i$xix) +i^(ln(a)x(x) 



2 4 7T 2 



(z^ 2 - ~) J $ - e 1 ^ / '' cfo/a'j [^(iO^Ar/ 2 ) + l] sm(kArj) 



7T1 



cos (kAr]) si(2/cAr/) + — + sin(A;Ar/) ci(2/cAry) — 7^ — In 



kAr] 



x(v') 



+ ln(a)i$ — e 

rv 



k 



ik-x 



drj'a' sm(kAr])x(r]') 



drfa! ln(a') sm.(kAr])x(f]') 
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/ 2 tfV 



Ak-x 



drj'a | ln(K 1 H 2 At] 2 ) sm(kAr]) - kAr/ cos(kAr]) 
+ (kAr ] ) 2 akAr ] )}x(r ] ') 



+ m(a)i^ p-e 



lh-x 



drj'a' 2 sm(kArj) — kArj cos(kArj) x(jl') 



1 



+i$ tW*'*/ dr/'a w ln(a') sin^Af?) - kAr] cos(/cAr/) x (t/) !> = , (28) 



fc 3 



/2 



where we made use of (writing z = kAr]) 

(d 2 + £; 2 )z 2 £(z 

(<9 2 + fe 2 ) 2 ln(az) ( sin(z) — z cos( 



2fc 2 {-cos(z) si(2z) + ^ + sin(z) ei(2z)-7 E -ln (|) 



d sin(z) — z cos(z) 



2fc 2 {[2 ln( 



az) + 1 



sin (2) + 



dz z 
d sin(z) — z cos(z) 



dz 



} • (29) 



The first term in (|28|) is the conformal vacuum contribution and thus does not grow during inflation. 
Upon conformal rescaling, this part of the self energy can be written as a function of x — x' , such 
that it can be analysed in momentum space, in which it appears local, as it is standardly done. The 
second contribution in (|28[) is the conformal anomaly, which is local and hence cannot be further 
simplified. The anomaly grows logarithmically with the scale factor, and hence it is subdominant 
when compared with the third term, which grows as aa! and aa']n(aa'). This contribution is the 
dominant one in the limit when s = (m 2 + ^R)/H 2 — ► 0, since it contains terms that go as 1/s, 
while the last (fourth) term is of order s°, and we shall not analyse it further. 

At this stage, we note that the fermion apparently does not experience the effects of a de Sitter 
invariant thermal scalar bath at a temperature Th = H/2n. Since this thermal bath should also 
be present for a conformally coupled scalar its interaction with the fermion has to be described 
either by the conformal vacuum or the conformal anomaly contribution. While the former term 
does not involve H and therefore contains no information about de Sitter expansion, also the 
anomaly cannot mediate such a phenomenon due to its manifest breaking of de Sitter invariance. 

Keeping - in accordance with the above discussion - only the first three terms, Eq. (|28|) simplifies 

to 



(i7°<9 -J-k)x(k,r]) + ln(a)(i7°<9 - 7 ■ k)x(k, 77) + (i7°<9 -7 • fc)(ln(a)x(fc,77)) 



f 2 H 2 a ( 2 1 



2 4 



7T 



1 



- ) { (i7 U <9 - 7 • k) - J dr/o'l [ ln(KH Ar] ) + lj sm(kAr]) 

- cos(kAri) si(2kAr]) + ^ + sin(fcAr?) ci(2kAr]) - 7^ - In (pp) 



+ In(a)(i7°flb-7-fc) 



rt 1 



dr] a sin(kAr])x(r]') 
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+ (ij do — 7 • k) — [ drf a 1 \n(d)sm.{k£\ri)x(r}') 



0. 



(30) 



Before we proceed, we note that 



(i7 (9 - 7 • k)r 0(Ar]) sin(A;Ar/) 



(31) 



is the retarded Green function of the operator (i7°(?o — 7 • such that the leading order term 
(contaning the ln(K)) 0(1/ s) satisfies a local second order differential equation. The local equation 
is obtained by acting with the operator — a(i7°9o — 7 ' k)/a upon Eq. (jHH)) (of course, the generic 
terms of the order O(s ) remain thereby nonlocal), 



(ad \ + k 2 + ia(a i) 7 5 ^)x(A:,??)+ ^ (a5 ^a + ln(a)fc 2 + ia(<9, 



ln(a) 



7 ZiA; x(fc,»7) 



+ ad -do + ^ + ialflb- )7 Ma)x(fc, ry) 



f 2 H 2 a 2 f 2 1\ f r /3 \ /3 \ i ri . , . , , 



+ -Sb / dT]'a'ln(H 2 A7] 2 )x(v') + - drj'a 



2 n ,cos(kAr]) - 1 



Ar] 



+ -(doln(a)J J drj'a' cos(kAt])x(r]') 

+ - (d H a i) ih-f P drfa' sin(A;Ar/)x(?? / ) | = ' 



(32) 



where we assumed that x 1S a helicity eigenvector, 7°7 • kx(k,rj) = khx(k,rj). 

From this form, we infer that the conformal anomaly combines neatly with the tree level operator 
to yield 



fZ fZ q 

(l + Ha)) (do + k2 ) -oH(} + [Ha) " o]) a ° " ' iaHhk { 1 + 



16vr 2 



ln(a) 



(33) 



such that the conformal anomaly correction to the fermionic wave function may be neglected when 



ln(a) < 



16tH 



p > 



(34) 



while from ()33j) . one sees that a different definition of the scale factor a can always be absorbed in 
the field strength renormalisation. 

The leading order contribution 0(s^ 1 ) to the modified Dirac equation Q32JI can be easily ex- 
tracted from Eq. ()16l) . Dropping the anomaly term we get 



(ad -d + k 2 + ia(d -)hk^ 5 )x(k,ri) + ^ 



3f 2 H 2 


r h 2 


8tt 2 


m 2 +£R 



X(k,v) =0. 



(35) 
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This is to be compared with the corresponding expression for a free massive conformally rescaled 
Dirac-fermion in an expanding FLRW space, a special case of which is de Sitter inflation, 

1 d (X d el 

—a(\d—am)-{\(/)+am)x{x) = d x(v)+k x(v) — —d v x( r ))~^-^—h\k.\'j 5 x( r l)+ a2rn2 x( r )) = 0- (36) 
a 1 a a 

A comparison of Eqs. (|36|) and (|35|) suggests that the one-loop radiative corrections due to a Yukawa 
coupling to a light scalar field generate an effective mass for the fermion, 



3f 2 H* 
^ ~ 87r 2 (m 2 +£R) 



,2 _ U J ±x - ,.2tt2 



H H . (37) 



This is the main finding of this paper. A comparison with the corresponding result for the mass 
niy of a photon P, 0] coupled to a light charged scalar (e denotes electric charge), 

2 3e 2 H A 

suggests that the dynamical mass generation is a generic feature of de Sitter inflation. However, 
as we argue in sections |1] and El below, a more subtle analysis is required to fully apprehend all of 
the ramifications of the fermionic mass generation mechanism. 

In the original work Ref. an instability of the fermionic wave function was claimed, which is 
due to an incorrect sign in front of the fermion self energy, and hence unphysical. When the sign 
error is corrected |l(, the wave function exibits a behaviour that is similar to the one reported here 
and which is consistent with mass generation in the case of a Yukawa coupling to a massless scalar 
field. 

4. INFLATIONARY FERMION DYNAMICS 

In this section, we derive analytic solutions to Eq. (|3"5)) and discuss their properties as well as 
some implications. We first note now that equations (|35H36j) can be reduced a first order system 

e=d/d V ), 

\L h — hkLh = am^pRh , 

\R' h + hkRh = am^Lh , (39) 

where we have decomposed the 4-spinor x into a direct product of chirality and helicity 2-spinors 
as follows, 




X =\ ® & , ( 4 °) 
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and where £/, is the helicity 2-eigenspinor, 

kh = k-&^ h = K h . (41) 

It is useful to consider the spinors 

u±h = ^= , (42) 

sucht that Eqs. ()39|) turn into 

\u' +h - am^u +h = hku_ h , 

iu'_ h + am^u-h = hku + h ■ (43) 
These can decoupled into two second order equations, 

u± h + (k 2 + a 2 m 2 b )u± h ± i{am^)'u± h = , (44) 
which in de Sitter space read 

/ I _ (I ii^ 2 \ 
u'i h + k 2 + 4 V2 , HJ )u ±h = 0. (45) 



This is Bessel's equation for the index 



2 1 H 

and has the normalisable solution 



ry 2 



1 fAC\ 

^±=0^1^ ( 46 ) 



u +h = 



-^< ) (-^)- (47) 



Making use of Eq. I|43|) and the recursion relation 



H^z) = -£-2#>(*) + -H®(z) (i = 1,2), (48) 



we obtain the second solution, 



U- h = hJ*"-J-^HM(-kn), (49) 



such that both solution are normalised to one, 



£h ±/l | 2 = i. (50) 
± 



This normalisation condition follows from the requirement that (in spatially homogeneous space- 
times) the vector current for each helicity sector must be conserved (see e.g. Ref. [11]), d^j^ = 0, 
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a consequence of Noether's theorem applied to phase transformations of the fermion field. It 
therefore holds for any fermion with time-dependent mass term, and its physical meaning is charge 
conservation. The analogous condition for a scalar field is the conservation of the Wronskian. 



We therefore introduce the phase space charge density associated with the vector current (j° 



Jd 3 kf 0h /(27rf), 



foh = \Lh\ 2 + \Rh\ 2 = \u+h\ 2 + \u-h\ 2 = 1 ■ 



(51) 



which is normalised to the vacuum fluctuations contribution. In order to check explicitly that our 
solutions satisfy Eq. (|50|) . we shall make use of (f4*8|) and the following identities 12] 2 



1 



Vj, = V- 



H®,(z) = J™H®(z), (i = l,2) 



(2), 



Note first that Eqs. (|5l|). (gTJ) and (@1J) imply, 



(52) 



foh 



it z r 
Ik 



iW*){v+-»-)H${z)HW{z) + e- i ^ /2){u +- v - ) H^}{z)H^}{z) \ , (z = -krj > 0) . (53) 



Making use of Eqs. ()52|) . we then obtain 



foh 



. TTZ 



Upon making use of the recursion relation (|48[) and the Wronskian, 



H^(z) = ^H^{z) + -H®(z) (i = 1,2) 



d 



dz 



this immediately reduces to 



(54) 



(55) 



foh = 1 > 



(56) 



proving thus Eq. (f5T|) . 

Even though there is no charge generation during de Sitter inflation by the fermion mass gener- 
ation mechanism under consideration, there is particle number generation, which can be calculated 
following Ref. 



1J as 



ii 



kh 



(57) 



2 The last identity in 152H is incorrectly stated in Ref. |l2(|. 
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where u>k = ^k 2 + a 2 m 2 b and f\h, fih and fsh are the scalar, pseudoscalar and pseudovector 
densities on phase space, respectively, defined as 

f lh = -2U[L h R* h ] = -(\u +h \ 2 - \u. h \ 2 ) , (58) 
hh = 2%[L* h R h ] = 2Q[u +h u*_ h ] , (59) 
hh = \Rh\ 2 - \L h \ 2 = -2M[u +h u*_ h ] . (60) 

It is known 11] that the particle number density on phase space n^ h lies between and 1, in 
accordance with Pauli blocking. Note however, that for adiabatic particle production, which is the 
case at hand, one should interpret n^ h with care, as it may have very different properties than a 
corresponding particle distribution in the time-independent case. For example, an Unruh detector 
coupled to a scalar field in an expanding background is apparently insensitive to the corresponding 
scalar expression for n? h 

on . As n% h is derived by diagonalisation of the Hamiltonian, a quantity 
with a clear physical meaning is the energy per mode, Qr, = u^n^, which defines the total energy 
density (the vacuum contribution is subtracted), p = J2h=±l I d 3 kQ,^ h /(27r) 3 . 
For the situation at study, Eq. ()57l) reduces during de Sitter inflation to 
1 1 



n 



kh 



2^1 + ^/z- 



-.[2h^[u +h u*_ h ] + -{\u +h \ 2 -\u^ h \ 2 )) 



(61) 



which we now discuss in the super- and sub-Hubble limits, respectively. Based on the small 
argument expansion of the mode functions, 



U+h 



U-h 



2tt 
\h 



e^ + F(, + ) It 



e^I>_)(^ 



-v+ 



1 + e~ iw + 
1 + e~ i7TU - 



r(-i/+) (z\2v+ 
r(-i/_ 



-)' 

2) 
2< 



it is easy to obtain the late time limit, z = 

kri^O- 1 



n 



r(i/_) 

-kr) — » of expression (|6*Tj) . 
- tanh(7T/u) 



l + 0{z 2 ) 
l + 0(z 2 ) 



(62) 



kh 



1 



+ 0(zlm») 
+ 0(zln(z)), 



(63) 



which is reminiscent of a Fermi-Dirac distribution of (nonrelativistic) particles with (de Sitter) 
temperature, Tjj = H/2ir and an energy, E 1 ^ ~ m^. For light particles, <C H, we have, 
nr. ~ 1/2, while for heavy particles, S> i7, the population density is - as expected based on e.g. 
adiabatic analysis - exponentially suppressed, n^ h ~ exp(— 2-Km^/H). We have thus established 
a mass generation mechanism during inflation, which is responsible for particle production on 
super-Hubble scales, 

k 



(64) 
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according to the distribution (|63[). 

It would be incorrect to infer from the distribution (jHSJ) that the population density is thermal, 
since that distribution applies only for the infrared modes, which satisfy Eq. (|64jl. To get more 
information about the spectrum, we need to investigate the momentum distribution, which can 
be obtained in the ultraviolet limit z 3> m^/H by an expansion of the particle number (|61[) in 
my /(Hz), which turns out to be 

,2 



rn 



ii 



kh 



16z 4 H 2 
m%H 2 
(2fc/a) 4 



+ o 



+ o 



H 4 z Q J 
Jk/af 



(65) 



Hence fermions do not acquire an exponentially falling thermal distribution in the momentum k, 
as one would expect if they were interacting with a thermal bath of scalar particles. 
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FIG. 1: Fermion particle number n^h as a function of the fermion mass, m^p/H. 4 snapshots are shown, 
k/(Ha) = 1 (first Hubble crossing), k/(Ha) = 0.1,0.01,0 (end of inflation). 



In FIG. Q we show the evolution of the particle number, n^ h = n^ h (m^/H), during de Sitter 
inflation. Four snapshots are shown: k/(Ha) = 1,0.1,0.01,0. Note that most of the particles are 
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created after the first Hubble crossing at k/a\ x = H. At the end of inflation (k/(Ha) ~ 0) the 
particle number density approaches the one in Eq. (|63|) . 

An interesting question is what happens to the particles produced during inflation in the subse- 
quent epochs of preheating, radiation and matter domination, when we expect the fermionic mass 
to dissolve. This and related questions are the subject of a forthcoming publication. 



5. DISCUSSION 



As main result of this paper, we have presented the effective fermion mass (jl~>7j) . dynamically 
generated in de Sitter space. From the thermal features of de Sitter space, and especially from the 
periodicity of the scalar Green function in imaginary time , one might expect that additional 

degrees of freedom coupled to the scalar field thermalise. This would be in analogy with an Unruh 
detector, which has a thermal response function and, as a consequence of the principle of detailed 
balance, equilibrates at de Sitter temperature. However, we find no indication of scatterings from a 
thermal bath or even thermalisation experienced by the fermion through interaction with the scalar 
field. The apparent reason is, that for an Unruh detector, it is assumed that its internal dynamics 
is governed by its proper time and that detailed balance holds for the detector [141 ] . One might 
question whether this requirement may in principle be realised by an experimental device without 
disturbing the curved background to an extent which spoils the measurement. This is a concern 
which does not apply to the field theoretical investigation conducted in this work. Nonetheless, 
the mass generation mechanism fits into the Unruh detector picture when relating it to the Lamb 
shift of the detector's energy levels [3], which can also be interpreted as a radiative correction. 

Finally, let us discuss in more detail, how the dynamical mass term is related to the familiar 
Dirac fermion masses. We first note that we can describe the mass generation effect by a nonlocal 
effective action, 

,2 -, 

+ O(s ) + S anoma i y , (66) 



S off = d x 



where denotes the fermion mass given in Eq. ()37|) . l/(i<^) is the (retarded) nonlocal operator 
whose kernel is given in (|31|) and 

f 2 



'anomaly 



64vr 2 



a 3 / 2 ip ( ln(a)i$ + $ ln(a)) a 3 / V 



(67) 



is the effective action for the anomaly. 
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The question is then, what is the difference between the dynamics induced by an ordinary 
fermion mass term, as implied by the standard action, 

Smassivo = / d 4 x a 3/2 ipi^a 3/2 ip - a 4 m^ip , (68) 

and the dynamics generated by the nonlocal effective action Q66JI . As we have shown above, these 
two actions lead to identical second order evolution equations for the chiral densities, Lh and Rh 
(and likewise for their linear combinations, u±h)- Yet there is a difference: the nonlocal action (|66|) 
conserves chirality, while its local counterpart violates chirality. According to the local action (|68[). 
the left and right handed densities evolve according to the first order system (|39|). whose graphical 
representation is shown in terms of mass insertions in FIG. [2j On the other hand, the nonlocal 
action 1)66(1 cannot flip chirality, and its proper graphical representation is shown in FIG. |21 This 
is the sense in which the fermion mass 1)37)1 is not the genuine Dirac mass, which couples left and 
right handed amplitudes. One may think of the equation of motion for the L-handed (R-handed) 
fermions as being governed by the system of equations ()39)) . but one should keep in mind that 
the fermions of opposite chirality are "ghost" fermions and do not comprise physically measurable 
states. This is to be contrasted with the photons in de Sitter background, which acquire a mass 
through interactions with a scalar medium, and thus an additional longitudinal physical degree of 
freedom Q,[J. 

L R L R 
X X X 

Hfy Hfy 1% 

FIG. 2: Mass insertions for the local fermion Dirac mass term 1)68)1 . These insertions violate chirality. 

L L L 
X X 



my 

FIG. 3: Quadratic fermionic mass insertions corresponding to the nonlocal effective action 1)66)1 . which 
preserve chirality. 

Putting these findings together suggests dynamical mass generation for interacting quantum 
fields as a phenomenon generically occuring in de Sitter space. Therefore, it is well conceivable 
that this mechanism may cause effects during cosmic inflation which lead to observable signatures. 
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